Periodic-orbit bifurcations as the origin of nuclear deformations 
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Semiclassical analysis of shell structures in realistic nuclear potentials are presented using periodic- 
orbit theory. We adopted r" potential model and examined classical-quantum correspondence us- 
ing Fourier transformation technique. Spin-orbit coupling is also taken into account in the model 
Hamiltonian. Gross shell structure for a certain combination of surface diffuseness and spin-orbit 
parameters are investigated and its relation to pseudospin symmetry is discussed. Analysis of su- 
perdeformed shell structure in realistic model is also presented. 
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1. SHELL STRUCTURE AND PERIODIC 
ORBITS 

Nuclear deformations are intimately related with shell 
structures in single-particle energy spectra of deformed 
Hamiltonian. Using the semiclassical theory, quantum 
level density g{E) can be represented in terms of classical 
periodic orbits, and one obtains the trace formula[l, 2] 
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go{E) is average part of the level density, and the os- 
cillating part is expressed as the sum over all periodic 
orbits f] in corresponding classical Hamiltonian system. 
3/3 = §^p-dr is action integral along the orbit /3, and Vfj 
is Maslov phase determined by the number of conjugate 
points along the orbit. 

In the trace formula, each periodic orbit contribution 
is an oscillatory function of energy since action Sp{E) 
is an increasing function of E. The energy scale of this 
oscillation is given by 6E « 2Tih/Tp. From this relation, 
one can see that short periodic orbits (having small T^) 
is associated with gross structure (having large 5E). 

The amplitude factor Afj in the trace formula (1) has 
significant dependence on the stability of the orbit. In 
standard stationary-phase approximation, it is propor- 
tional to the stability factor; 
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where Mp represents symmetry-reduced monodromy 
matrix[3], which characterizes linear stability of the peri- 
odic orbit. Varying external parameters in the Hamilto- 
nian, each periodic orbit continuously changes its prop- 
erties, and one of the eigenvalues of Mg might eventu- 
ally approaches 1. At this point, a continuous family of 
quasi-periodic orbits appear in neighborhood of the peri- 
odic orbit (3 in direction to eigenvector belonging to the 
above unit eigenvalue of Mp. It usually accompany the 
appearance of new periodic orbit from that local fam- 
ily (or, inversely, disappearance of another periodic or- 
bit into the family), namely bifurcation of periodic or- 
bit occur at this point. One should also note that this 



periodic-orbit bifurcation is related with the restoration 
of local dynamical symmetry in neighborhood of the peri- 
odic orbit. Since these quasi-periodic orbits make coher- 
ent contribution in the periodic orbit sum, we can expect 
a significant enhancement of shell effect in quantum level 
density. Actually, our previous works [4, 5] clearly show 
that the bifurcations of short periodic orbits play essen- 
tial role in emergence of gross shell structures. The di- 
vergence of stability factor (2) at bifurcation point is due 
to the breakdown of standard stationary-phase method, 
and can be remedied by using suitable approximation for 
trace integral such as uniform approximations[6, 7] and 
improved stationary phase method [8]. 



2. THE r" POTENTIAL MODEL 

For nuclei and metallic clusters, central part of the 
mean field is well described by Woods-Saxon potential. 
For stable nuclei, Woods-Saxon potential can be further 
approximated by much simpler potential with r" radial 
dependence [9] 
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Here. 



1 + exp((r - R)/a) 
the parameter a controls the surface diffuseness. 



With typical value of Woods-Saxon parameters for nu- 
clei, we found that the quantum energy spectra of r" 
model nicely reproduce those of Woods-Saxon model up 
to about -10 MeV. 

Shifting the zero-point energy, we write the r" model 
Hamiltonian as 
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The advantage of the r" model is its scaling property. 
Since equations of motion (EOM) are invariant under 
scale transformation 
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one has the same set of periodic orbits in all energy sur- 
faces. The action integral has simple energy dependence 
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where we define scaled energy £ as 
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We also define scaled period by Tp = Sf3{E)/h£, which is 
a scale-invariant parameter proper to the periodic orbit 
f3. Then we obtain semiclassical level density with scaled 
energy 

g{£) = g{E)^ = go{£) + M^)^os {£Tfs - up). (8) 

This simple form is convenient for Fourier analysis. The 
Fourier transform of the trace formula 
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is a function having peaks at scaled period of classi- 
cal periodic orbits. It suggests that wc can extract in- 
formation on classical periodic orbits from the Fourier 
transform of quantum level density 
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where £i is related with quantum energy level Ei via 
Eq. (7). 

For description of nuclear shell structure, it is very 
important to incorporate strong spin-orbit coupling in 
the mean field. One can also apply periodic orbit the- 
ory to a system with spin-orbit coupling, if we uti- 
lize SU(2) coherent state path integral method[10 12]. 
It introduces classical variable corresponding to a spin 
degree of freedom. Let us write spin vector in polar 
coordinate as s = s(sin cos <p, sin i? sine/?, cos with 
s = h/2. One can choose spin canonical variable as 
{qs,Ps) = (95, s cost?) and obtains the Poisson bracket 
relations {si, Sj}p,B. = eijkSk, which exactly correspond 
to the commutation relation of quantum spin operators. 
The Hamiltonian with spin-orbit coupling is written as 
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In order to make use of scaling property of the central 
potential, we adopt spin-orbit potential 



(12) 



and choose the diffuseness parameter a' = a/2 + 1. 
With this choice of spin-orbit coupling, the classical EOM 
are invariant under transformation (5) for some special 
classes of periodic orbits whose spin does not change its 
orientation; s — 0. For axially symmetric or reflection- 
symmetric shapes, such orbits exist as exact solutions of 
the EOM. Using the Hamiltonian (11) with suitable a 
and K, we can reproduce realistic shell structures both 
for spherical and deformed nuclei. 




FIG. 1: Scaled periods Tg of classical periodic orbits plotted 

as functions of spin-orbit parameter k. Open circles repre- 
sent the bifurcation points. The background gray-scale image 
show the Fourier transform of quantum scalod-cnorgy level 
density; the Fourier amplitude |F(T;ft)| have peaks in dark 
regions. 
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FIG. 2: Singlc-particlo level diagram against spin-orbit pa- 
rameter K. Solid and broken curves represent even and odd 
parity levels, respectively. Open circles indicate the crossing 
points of pseudo spin-orbit partners. 
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FIG. 3: Oscillating part of coarse-grained quantum level density is compared with semiclassical trace formula. In the top panel, 
the shortest circular orbit and oval orbit (2,1) are considered in the periodic orbit sum. In the bottom panel, longer circular 

orbit and triangular orbit (3,1) arc also taken into account. (In this plot, amplitude Ap and Maslov phase in trace formula 
(8) are determined so that they best fit the quantum result.) The right panels display the classical periodic orbits, where the 
outermost circle represents the boundary of classically accessible region. 



3. NUMERICAL RESULTS 

In Fig. 1, we compare scaled period rg's of classical 
periodic orbits with Fourier transform of quantum level 
density. Absolute value of Fourier transform |-F'(t)| are 
displayed by gray-level image and it exhibits peaks at 
dark regions. We found that all important peaks of 
Fourier amplitude are located at scaled periods T/j of 
classical periodic orbits, displayed by curves, and they 
arc significantly enhanced around the bifurcation points. 
Especially, the Fourier amplitude corresponding to cir- 
cular and (3,1) orbits are strongly enhanced around the 
bifurcation point k ~ 0.05. It implies a significant growth 
of gross shell structure around k ~ 0.05. 

Figure 2 shows the single- particle level diagram against 
spin-orbit parameter k. As we expected from the Fourier 
analysis, remarkable shell structure emerge at about k ~ 
0.05, with large energy gaps corresponding to well known 
nuclear magic numbers 50, 82 and 126. One should also 
note that pseudo spin-orbit partners degenerate at al- 
most same values of k. These degeneracies are responsi- 
ble for the occurrence of the above mentioned large en- 
ergy gaps. This implies that pseudospin symmetry[13] 
is approximately restored in this region. As we stated 
in §1, local restoration of dynamical symmetry is usually 
accompanied by bifurcations of classical periodic orbits. 
Therefore, we can expect that the pseudospin symmetry 
might be connected with periodic orbit bifurcation. 

Pseudospin symmetry arise for a certain combination 



of spin-orbit and orbit-orbit force parameters in Nilsson 
Hamiltonian. In the Nilsson model, orbit-orbit force is 
introduced to describe sharp surface of potential well and 
therefore connected with surface diffuscness. In r" po- 
tential model, bifurcations of periodic orbits occur for 
certain combinations of diffuseness parameter a and spin- 
orbit coupling parameter k. Therefore, it should be quite 
natural to consider that bifurcation of (3,1) orbit which 
occurs at k ~ 0.05 is a semiclassical origin of gross shell 
structure related with pseudospin symmetry. 

In order to see the importance of (3,1) orbit for gross 
shell structure, we analyze the coarse-grained quantum 
level density. In the left panels of Fig. 3, oscillating part 
of quantum level density coarse-grained with a smoothing 
width AS = 0.3 is plotted with dotted curve. In the peri- 
odic orbit sum, contribution of long orbits are suppressed 
by coarse-graining. For AS = 0.3, orbits with scaled pe- 
riods Tfs > 2tt/AS ~ 20 arc suflicicntly suppressed. The 
solid curve in the upper panel of Fig. 3 represents the 
semiclassical level density where only shortest circular 
orbit (whose spin is anti-parallel to the orbital angular 
momentum) and oval orbit (2,1) (see the right upper part 
of Fig. 3) are taken into account. It apparently fail to re- 
produce the gross pattern of quantum result. Next we 
take into account the contribution of longer circular or- 
bit (whose spin is parallel to the orbital angular momen- 
tum) and triangular orbit (3,1) (see the right lower part 
of Fig. 3). Then we can nicely reproduce the gross struc- 
ture of quantum level density as shown in the lower panel 
of Fig. 3. Thus we see that this bifurcating orbit actu- 
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structure. 

Using the r" model with spin-orbit couphng, we also 
analyze the semiclassical origin of superdeformed shell 
structure. FigTirc 4 shows thc> Fourier transform of 
level density as function of deformation parameter 6 = 
3(77— l)/{2ri + 1), connected with axis ratio rj = Rz/R±. 
We found significant enhancement of Fourier amplitude 
at (5 = 0.5 ~ 0.6, corresponding to the bifurcations of 
periodic orbits shown in Figure 5. With increasing defor- 
mation, oval orbit lying around the short diameter occur 
period-doubling bifurcation and new orbits (4,2, l)a and 
(4,2, l)b emerge at S ~ 0.5, and then they submerge into 
the oval orbit lying around the long diameter at 6 > 0.6. 
This situation is quite similar to the case of the model 
without spin-orbit couplings [5]. 



FIG. 4: Scaled period of classical periodic orbits plotted as 
functions of deformation parameter S. Open circles represent 
the bifurcation points. In the background, Fourier transform 
of level density |F(t, is plotted by gray-scale image. The 
Fourier amplitude exhibits peaks at dark regions. 
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FIG. 5: Periodic orbits responsible for superdeformed shell 
structure in r" model with a = 3.0 and k = 0.05. 



ally play essential role in characterizing the gross shell 



4. SUMMARY 

We applied periodic orbit theory to a realistic nu- 
clear mean field potential model with spin-orbit cou- 
pling. The scaling property of r" potential model en- 
ables us to investigate quantum-classical correspondence 
using Fourier transformation technique. We found nice 
quantum-classical correspondence, and significant roles 
of periodic-orbit bifurcations are clarified, also for our 
systems with spin-orbit coupling. Especially, (3,1) orbit 
bifurcation which occur for a certain combination of dif- 
fusencss and spin-orbit parameters can be regarded as 
the semiclassical origin of gross shell structure related 
with pseudospin symmetry. Superdeformed shell struc- 
tures are also examined and we found significant roles of 
bifurcations of periodic orbits which are similar to those 
in the model with no spin-orbit couplings. It implies that 
the semiclassical origin of superdeformed shell structure 
is essentially accounted for by models without spin-orbit 
couplings. 
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